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OUTLINE  OF  METHODS  FOR  DESIGN  OF  SUPERCONDUCTING  TURBOGENERATORS 

Ludwik  Antal  and  Jacek  Uczkiewicz 

Most  Important  Notation 

Ag  =  specific  electric  loading  of  stator 
B  =  magnetic  induction 

=  dimensional  coefficients 

z,v,u,L 

D  =  diameter 
i,I  =  current 

j  =  current  density 
kw  =  filling  factor 
1  =  length 

lga  =  calculated  length  of  armature  winding 
n  =  rotational  velocity 
r,R  =  radius 

S  =  apparent  power 
T  =  time  constant 
Y  =  electric  conductivity 
A.  =  shield  thickness 
X  =  relative  permeance 
u  ,y*  °  magnetic  permeability 

o 

$  -  phase  shift  angle 
w  =  angular  frequency 


Introduction 


Work  on  synchronous  generators  with  a  superconducting  excitation  winding 
is  going  on  in  the  most  important  scientific-research  centers  in  many  coun¬ 
tries.  The  programs  of  activities  undertaken  in  this  field  are  being  realized 
very  persistently,  and  everything  indicates  that  superconducting  turbogenera¬ 
tors  already  occupy  a  permanent  place  among  electric  machines.  The  efforts 
of  scientific  teams  are  currently  being  concentrated  on  several  very  difficult 
problems.  The  latter  include  heat  calculations  of  a  rotor  comprising  a  super¬ 
conducting  winding  cooled  by  liquid  helium,  electromagnetic  calculations  of 
the  excitation  winding  shielding  system,  dynamic  stability  studies,  strength 
calculations  of  the  supporting  structure  of  the  rotor  and  the  structure  clamp¬ 
ing  the  slotless  armature  winding.  Concentration  of  effort  on  selected 
problems  is  taking  place  at  the  expense  of  a  general  design  method,  resulting 
in  the  loss  of  relationships  between  the  results  of  analyses  of  individual 
phenomena.  The  design  of  power  engineering  units  (.output  exceeding  1  GVA) 
necessitates  the  development  of  computational  methods  which  are  suitable  for 
an  analysis  of  phenomena  that  are  negligible  or  unknown  in  conventional 
machines.  On  the  one  hand,  these  methods  must  describe  with  the  greatest 
possible  accuracy  physical  phenomena,  and  on  the  other,  they  must  be  inter¬ 
related  to  the  maximum  possible  degree.  This  follows,  among  other  things, 
from  economic  prerequisites— the  cost  of  building  a  unit  with  output  on  the  order 
of  magnitude  1  GVA  is  too  high  to  permit  the  application  of  traditional  methods 
for  improving  the  design.  Therefore,  a  need  exists  for  formulating  an  algorithm 
which  will  determine  rationally  basic  design  magnitudes  of  superconducting  gene¬ 
rators  on  the  basis  of  a  synthesis  of  the  results  of  analyses  of  individual 
phenomena. 

Description  of  Algorithm 

The  first  actiyity  during  the  design  of  a  synchronous  superconducting  gene¬ 
rator  is  selection  of  the  structure.  Other  considerations  pertain  to  a  machine 
with  the  design  envisioned  for  this  type  of  converter  in  the  output  range  0.3  to 
3  GVA  I6J .  Fig.  1  shows  the  basic  elements  of  the  electromagnetic  system  and 
their  mutual  arrangement.  The  main  purpose  of  using  the  algorithm  is  the 
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determination  of  basic  geometric  dimensions  which  will  correspond,  as  accurately 
as  possible,  to  the  optimal  solution.  These  dimensions  are  determined  by  deter¬ 
mining  coefficients  defined  by  their  ratios.  The  following  dimensional  coeffi¬ 
cients  are  used  (.Fig.  l) 
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[0  -  calculated;  c  =  armature;  w  =  filling;  v  =  calculated; 
e  a  shield;  cz  =  frontal] 


Fig.  1.  Electromagnetic  system  of  superconducting 
generator 

1.  superconducting  excitation  winding 

2.  rotor  shielding  system 

3.  slotless  stator  winding 

4.  magnetic  shield  of  stator 
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!Hie  next  steps  in  formulating  the  design  algorithm  are:  determination  of 
the  structural  equation  of  the  machine  and  setting  up  criteria  for  selection 
of  dimensions.  Introduction  of  these  criteria  in  the  structural  equation  re¬ 
duces  the  number  of  successive  approximations  in  the  endeavor  to  design  a 
suitable  optimal  structure. 

In  the  preliminary  design  stage,  these  criteria  are  determined  by  assuming 
great  simplifications  in  detailed  computational  methods ,  which  allows  one  to 
narrow  down  the  range  of  possible  variations  of  some  dimensions.  The  next 
step  in  the  algorithm  is  proper  sequencing  of  detailed  electromagnetic,  mech¬ 
anical,  and  heat  computational  methods  and  their  assignment  to  a  common  base 
constituting  a  description  of  the  electromagnetic  fields  in  the  machine.  Qua¬ 
lity  measures  of  the  design  are  verified  on  the  basis  of  the  results  of  de¬ 
tailed  calculations  and  if  necessary  the  dimensional  coefficients  adopted  in 
the  initial  stage  are  corrected.  The  algorithm  for  this  procedure  is  presented 
in  Fig.  2  in  the  form  of  a  block  diagram. 

Structural  Equation  of  Machine 

The  fundamental  structural  equation  of  a  superconducting  machine  is  the 
same  as  that  for  a  conventional  machine  [UJ 

^.kX9B„A.L,':»%  (2) 
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The  equation  in  this  form  is  not  suitable  for  the  design  of  superconducting 
generators  for  two  reasons:  first,  it  does  not  take  into  account  all  basic  dimen¬ 
sions  of  the  machine  (.Fig*  1);  second,  the  discrepancies  resulting  from  the  sim¬ 
plifying  assumptions  made  during  its  derivation  are  so  great  for  the  machines 
under  consideration  that  the  equation  should  only  be  treated  as  a  preliminary 
descriptive  equation.  The  complete  structural  equation  was  obtained  by  relating 
the  power  of  the  superconducting  machine  to  its  basic  dimensions,  current  den¬ 
sities  of  windings,  and  angular  velocity  X 2 J : 
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Fig,  2.  Block  diagram  of  algorithm  determining  basic 
dimensions  of  synchronous  superconducting  machine 


a.  Ratings 

b.  Selection  of  structure 
and  materials 

c.  Criteria  for  selection 
of  dimensional  coeffi¬ 
cients/initial  calcula- 
tion  of  dimensions 

d.  Structural  equation 

e.  Electromagnetic  calcula¬ 
tions 

f.  Heat  calculations 


g.  Strength  calculations 

h.  Measures  of  quality  of 
design 

i .  Results 

J  .  Correction  of  dimensional 
coefficients 

k.  Yes 

l.  No 
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This  equation  can  be  simplified  by  assuming  some  quantities  to  be  constant 
for  example,  the  filling  factors  of  properly  made  windings  [2}  are:  k  =  0.25 


'  nV.S'a'. 


k  _  =  0.625,  while  the  ratio  of  the  calculated  length  of  the  machine  1„  to  the 
wf  Ox 

calculated  length  of  the  armature  winding  1  is  approximately  equal  to  O.J. 

Ua 

In  addition,  assuming  p  =  1,  Eq.  (3)  takes  on  the  form 


S*  =  20.203  •  10"*  (1  -  <JV  l*3lsRAJ~  X 


cos2  pN  —  sm 


Eq.  (.4)  contains,  "besides  the  ratings  S^,  cos^,  the  ratios  of  the  radii 

of  the  windings:  c,  d,  x,  the  current  densities  of  the  windings  and 

the  product  R  1  .  The  unit  permeances  X  and  X,  are  functions  of  the  coeffi- 

ao  oa  ad  d 

cients  c,  d,  y  II].  The  fundamental  structural  equation  of  the  machine  is  ob¬ 
tained  by  transforming  relation  (.4)  in  such  a  way  that  the  outer  diameter  of  the 
armature  winding  is  expressed  in  terms  of  the  dimensional  coefficients  (l) 
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Criteria  for  Selection  of  Dimensions 


The  variability  ranges  of  the  current  densities  jfN  and  in  Eq.  (5) 

are  constant.  The  armature  current  density  is  determined  by  the  adopted  cool- 
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ing  system  and  is  in  the  range  (2.5  to  10) -10  A/m  .  The  excitation  winding 
current  density  is  determined  by  the  critical  magnitude  of  the  current  in  the 
superconductor  that  is  used,  which  for  the  given  material  depends  on  the  tem¬ 
perature  and  the  external  magnetic  field  strength.  The  characteristics  j  = 

=  f(.T,H)  determined  for  short  specimens  can  only  be  used  for  preliminary  cal¬ 
culations,  After  the  design  of  the  excitation  winding  structure  and  selection 
of  the  catalogued  conductor ,  the  excitation  current  density  assumed  earlier  can 

be  corrected  using  the  catalogued  characteristic  I  =  f(T,H).  The  critical  cur- 

cr 

rent  density  for  conductors  produced  at  the  present  time  are  in  the  range 

£  O 

(200  to  400)*10  A/m  at  T  =  4.2°K  and  H  =  397.5  A/m.  Assuming  that  the  rated 


excitation  current  density  is  equal  to  two-thirds  of  the  critical  density,  we 

obtain  a  practical  range  for  the  magnitude  of  the  excitation  current  density 
c  ? 

(130  to  270) *10  A/m  . 

The  dimensional  coefficients  c,  d,  x,  y,  L  follow  from  the  criteria  assumed 
in  the  design.  In  the  case  of  a  superconducting  turbogenerator,  these  criteria 
may  be: 

1.  dynamic  stability 

2.  critical  velocities  of  rotor 

3-  mechanical  stresses  in  stator  and  shielding  system  in  a  short-circuit 
surge  state 

4.  maximum  deformations  of  rotor  in  transient  states 
5*  time  constant  of  shielding  system 
6.  free  vibration  of  stator  core. 

Criteria  for  selection  of  dimensions  follow  from  a  detailed  analysis  of  the 
given  quantity  (dynamic  stability),  mechanical  stresses  in  rotor,  etc.,  which 
are  reduced  to  the  admissible  ranges  of  the  values  of  the  dimensional  coeffi¬ 
cients  using  the  properties  of  the  materials  intended  for  use  and  the  mechanical., 
electromagnetic  and  thermal  properties  of  the  designed  structure.  Because  many 
analyzed  criterion  quantities  are  mutually  opposed,  an  attempt  should  be  made  to 
determine  a  sufficiently  large  number  of  them  already  in  the  first  design  stage. 
This  allows  one  to  limit  substantially  the  number  of  design  stages. 

Detailed  Calculations  and  Quality  Measures 

Detailed  calculations  include  electromagnetic,  heat,  and  strength  calculations 
carried  out  for  steady  operation,  asymmetry  of  load,  and  short-circuit  surge.  The 
purpose  of  electromagnetic  calculations  is  above  all  a  determination  of  the  dis¬ 
tribution  of  the  field  in  the  machine.  This  constitutes  a  basis  for  calculating 
losses  in  individual  structural  elements ,  electrodynamic  volume  forces ,  electro¬ 
magnetic  parameters,  etc,  These  calculations  must  take  into  account,  as  needed, 
the  finite  axial  dimensions  of  machine  elements  and  the  effect  of  coil  outhangs. 


Heat  calculations  lead  to  a  determination  of  the  distribution  of  tempera¬ 
ture  fields,  heat  fluxes,  and  consumption  of  coolant  in  cooling  systems.  They 
also  allow  one  to  calculate  expansion  (dilatation)  forces  of  elements  operating 
at  great  temperature  gradients. 

Strength  calculations  allow  one  to  determine  mechanical  stresses  and  strains 
in  structural  elements  in  abnormal  states,  in  particular,  in  a  short-circuit 
surge  state,  while  rotor  strength  calculations  allow  one  to  improve  the  accuracy 
of  the  dimensional  criteria  already  in  the  initial  design  stage.  Detailed  cal¬ 
culations  should  in  fact  lead  to  relationships  between  parameters  of  the  inves¬ 
tigated  phenomena  and  the  dimensional  coefficients  of  the  machine.  However, 
because  of  the  complexity  of  these  phenomena,  in  many  cases  satisfactory  ana¬ 
lytical  solutions  have  not  yet  been  obtained.  In  particular,  this  applies  to 
the  distribution  of  fields  in  frontal  zones ,  stator  and  shielding  system  vibra¬ 
tions,  etc.  Until  these  phenomena  are  fully  described,  simplified  computational 
methods  must  be  used. 

The  quality  measures  of  the  design  include  two  groups  of  indicators.  Coef¬ 
ficients  of  strength,  heat  given  up,  etc.  constitute  the  first  group.  Verifica¬ 
tion  of  the  values  of  these  coefficients  and  their  comparison  with  the  admiss¬ 
ible  values  decides  whether  correction  of  the  dimensional  coefficients  is  neces¬ 
sary.  Operational  indicators  of  the  generator  (mass,  efficiency,  idle  run  char- 
aracteristic,  reactance,  etc.)  constitute  the  second  group  of  measures.  The 
values  of  these  parameters  constitute  a  basis  for  an  evaluation  of  the  design 
assumptions  and  possible  corrections  of  dimensional  coefficients  within  the 
ranges  determined  by  appropriate  criteria. 

Examples  of  Calculations 

To  illustrate  proper  selection  of  individual  dimensions  on  the  basis  of  the 
described  criteria,  we  selected  from  a  large  set  of  problems  the  problem  of  the 
critical  velocity  of  the  rotor  and  the  problem  of  the  time  constant  of  the  shield¬ 
ing  system. 
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In  regard  to  the  rotational  velocity  of  the  rotor  of  a  superconducting  gene¬ 
rator,  it  is  required  that  it  be  smaller  than  the  first  critical  velocity  or  that 
it  lie  between  the  first  and  second  critical  velocity.  Using  the  first  condition 
one  can  assume  [2j 

2 

njv  <  "*t 


Simplifying  the  problem  to  the  case  of  a  two-stage  rotor  made  from  nonmag¬ 
netic  steel  (Fig.  3),  we  obtain  a  relation  for  the  first  critical  velocity  as  a 
function  of  the  dimensions  of  the  rotor  [2j: 


nkt  =  3.04  •  10!  — - 

*C 


Substituting  condition  (6)  in  (7),  we  obtain  the  relation 

1 .  8.3 

Bm  (8) 

which  is  graphically  illustrated  by  the  plot  in  Fig.  3.  The  calculated  length 

of  the  armature  winding  ln  =  l.U*l  is  approximately  equal  to  the  length 

ua  8.  jqg an 

of  the  armature  sleeve,  whereas  the  outer  radius  of  the  armature  winding  R  . 

ai 

is  greater  than  the  outer  radius  of  the  rotor  (by  an  amount  equal  to  the  size 
of  the  gap  and  the  thickness  of  the  insulation  of  the  base  of  the  armature 
winding).  Assuming  R&^  =  1.2Rsq,  we  obtain  an  upper  bound  for  the  dimensional 
coefficient  L,  which  occurs  in  structural  equation  (5).  Hence,  during  selection 
of  the  value  of  L,  advantage  can  be  taken  of  the  graph  in  Fig.  3  and  of  the  in¬ 
equality 

L  <  0,833  d  — — 


In  regard  to  the  time  constant  of  the  shielding  system,  it  is  required  that 
it  be  1  to  2  times  greater  than  the  reciprocal  of  the  angular  frequency  of  the 
damped  field  13],  In  the  case  of  fluctuations,  this  is  equivalent  to  the  re¬ 
quirement 

T.>  0,2#  (m) 


Fig.  3,  Effect  of  rotor  dimensions  on  its  first 
critical  velocity 


1  /u  =  0.8;  R  /R  =  k-,  R 
r  c  so 

n^  =  3000  rpm 


front 


=  2; 


Key;  a'  Bfroat/Ro  =  2 
b .  rpm 


Therefore  the  damping  system  must  consist  of  two  shields :  a  damping  shield 
operating  in  the  surrounding  air  temperature  and  suppressing  fields  with  high 
frequencies ,  and  a  cryogenic  shield  operating  at  helium  temperature  and  ensuring 
an  adequately  large  time  constant  of  the  shielding  system.  The  system  of  shields 
must  be  protected  mechanically  by  suitable  nonmagnetic  steel  layers,  which  to¬ 
gether  form  a  multilayer  electromagnetic  system.  Determination  of  the  thickness 
and  subsequently  also  of  the  radii  of  individual  layers  is  of  fundamental  impor¬ 
tance  for  the  rotor  structure. 

The  time  constant  of  a  cylindrical  electromagnetic  shield,  without  taking 
into  account  the  yoke  of  the  stator,  is  calculated  from  the  relation 

T„  =  4  W* >Ri 


or  from  the  relation  taking  into  account  the  yoke  of  the  stator 


r.  =  r„ 


(12 


Fig.  U.  Time  constant  of  shielding  system  versus 

thickness  of  damping  shield  for  f  =  2  Hz,  p  =  1, 

V  =  1,  R  =  Q,4j4  m 
so 

1.  taking  into  account  stator  yoke 

2.  not  taking  into  account  stator  yoke 

3.  according  to  relation  Cll) 

4.  according  to  relation  (.12) 


Relations  (.11)  and  (12)  were  obtained  assuming  a  uniform  distribution  of 
the  density  of  eddy  currents  in  the  shield  and  for  f  -»*  0  1 5,].  In  reality,  the 
current  distribution  deviates  considerably  from  a  uniform  distribution  and,  in 
addition,  depends  on  the  frequency  of  the  field,  resistivity,  and  dimensions 
of  the  shield.  Therefore,  in  detailed  calculations,  the  time  constant  is  deter¬ 
mined  on  the  basis  of  the  actual  distribution  of  the  field  in  the  multilayer 
shielding  system.  We  considered  a  four-layer  system  consisting  of  an  external 


damping  shield,  a  supporting  structure,  a  vacuum,  and  a  cryogenic  shield.  The 
time  constant  of  the  system  is  determined  from  the  relation  [3.1 


T 

e 


i  r ,  . 

— arctan  AtB)  t,  /vT 


(13) 


in  which  T, (E  ),  Y,  (R  )  are  the  values  of  the  eigenfunctions  (r)  and  Y, (r ) 

1  so  1  so  11 

describing  the  distribution  of  the  field  in  a  cylindrical  multilayer  shield  [7]. 
Fig.  1*  shows  the  time  constant  of  the  shielding  system  as  a  function  of  the 
thickness  of  the  external  damping  shield. 


The  minimum  thickness  of  the  damping  shield  at  which  the  entire  system  has 
an  adequate  time  constant  is  obtained  from  the  plot  of  the  curves  presented  in 
Fig.  3.  This  -value,  combined  with  strength  calculations,  constitutes  one  cri¬ 
terion  for  selecting  the  outer  diameter  of  the  rotor. 


Conclusions 

The  idea  of  an  algorithm  for  the  design  of  basic  dimensions  of  electromech¬ 
anical  converters  applies  to  all  types  of  machines  in  which  problems  occur  that 
have  not  been  encountered  so  far  in  design  practice  and  for  which  a  structural 
equation  can  be  set  up  relating  the  ratings  to  these  dimensions.  The  design  of 
a  rational  structure,  that  is,  a  structure  substantiated  by  the  present  level  of 
knowledge,  is  the  reason  for  using  such  an  algorithm. 

The  degree  of  perfection  of  the  superconducting  synchronous  generator  de¬ 
sign  depends  on  the  number  and  importance  of  criteria  used  for  selecting  dimen¬ 
sions,  An  increase  in  the  number  of  criteria,  contributing  to  greater  computa¬ 
tional  accuracy,  is  limited  by  present-day  possibilities  of  an  analytical  descrip 
tion  of  phenomena.  The  selection  of  criteria  according  to  their  importance  en¬ 
tails  elimination  of  criteria  providing  constraints  on  dimensional  coefficients 
which  lie  within  the  ranges  determined  by  other  criteria. 


LITERATURE 


(11  Antal  L.:  Parametry  indukcylne  bezilobkowej  maszy- 
ny  synchronicznej.  Archtwum  Rleklrottchnxlei  (w  druku). 

(1|  A  n  t  a  1  L.,  Uctkitwicz  3. :  Metodyki  obliczert 

eieplnych,  elektrycznych  1  raechanicznych  dla  potrzeb 
komtrukcyjnych  prototypow  kriomaszyn  clektrycznych. 
Cx.  1,  Kaport  I-29/SPR27  Politechnlki  WrorlawskieJ 
Wroclaw  1919. 

I>1  Einstein  T.  H.,  Keirn  T.  A.:  Effect  of  field  winding 
electromagnetic  shield  of  the  excitation  requirements 
of  a  superconducting  alternator.  IEEE  Power  Eng.  Soc. 
Conf.  Pap.  Summer  Meet.,  San  Francisco,  Calif.  1975  A7 
59  79.5/1  —  At  56  79.5/9. 

N|  Miller  T.  J.  E.,  Hughes  A.:  Comparative  des  gn  and 
performance  analisys  of  air-cored  and  iron-core  i  syn¬ 
chronous  machines.  Proc.  IRE  vol.  24  February  19T7  No  2. 

(S|  Miller  T.  J.  E.,  Lawrenson  P.  J.:  Penetration  of 
transient  magnetic  fields  through  conducting  structures 
with  particular  reference  to  superconducting  a.c.  machi¬ 
nes.  Pro c.  IRE  Vol.  23  May  197S  No  5. 

(9)  Parker  3.  M.,  T  o  w  n  e  R.  A.:  Progress  In  supercon¬ 
ducting  turbine  generators  designe.  IEEE  PES  Winter 
Meeting  New  York  1977  A-77037-5. 

|7|  Ocikiewlcz  J.;  Cylindryczne  ekrany  wielowarstwo- 
we.  Archiwum  Elektrozechnitci  T.  XXVIII  :  1975  nr  4/110. 


